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Summary: We review a recently developed method, based on an exact 
auxiliary boson representation, to describe both Fermi liquid and non- 
Fermi liquid behavior in quantum impurity systems. Coherent spin and 
charge fluctuation processes are taken into account in a self-consistent 
way and are shown to include all leading and subleading infrared singu- 
larities at any given order of the self-consistent loop expansion of the free 
energy. As a consequence, for the SU(N) x SU(M) Anderson impurity mod- 
els the correct temperature dependence of the susceptibility is recovered 
over the entire temperature range, including Fermi liquid or non-Fermi 
liquid behavior below the Kondo temperature Tk- As a standard diagram 
technique the presented method has the potential to be generalized to 
correlated electron systems on a lattice. 



1 Introduction 



Highly correlated electron systems are characterized by a strong repulsion be- 
tween electrons on the same lattice site, effectively restricting the dynamics to 
the Fock subspace of states without double occupancy of sites. The prototype 
model for such systems is the Anderson impurity model, which consists of an 
electron in a localized level Ed < (called d-level in the following) with on-site 
repulsion U, hybridizing via a transition matrix element V with one or several 
degenerate conduction electron bands or channels jjj . Depending on the number 
of channels M, the model exhibits the single- or the multi-channel Kondo effect, 
where at temperatures T below the Kondo temperature Tk the local electron 
spin is screened (M — 1) or overscreened (M > 2) by the conduction electrons, 
leading to Fermi liquid (FL) or to non-Fermi liquid (NFL) behavior Q with 
characteristic low-temperature singularities, respectively. 

As perhaps the simplest model to investigate the salient features of correlations 
induced by short-range repulsion, the Anderson model plays a central role for 
the description of strongly correlated electron systems: In the limit of large 
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spatial dimensions g] strongly correlated lattice systems reduce in general to a 
single Anderson impurity hybridizing with a continuum of conduction electron 
states whose properties are determined from a self-consistency condition imposed 
by the translational invariance of the system [0. Quantum impurity models 
have received further interest due to their relevance for mesoscopic systems like 
single electron transistors or quantum point contacts || . Nonlinear conductance 
anomalies observed in the latter systems || have provided one of the strongest 
cases for the physical realization of the two-channel Kondo effect generated by 
two-level systems with electron assisted tunneling. 

The above-mentioned systems call for the development of accurate and flexible 
theoretical methods, applicable to situations where exact solution methods are 
not available. We here present a general, well-controlled auxiliary boson tech- 
nique which correctly describes the FL as well as the NFL case of the generalized 
SU(N)xSU(M) Anderson impurity model. As a standard diagram technique it 
has the potential to be generalized for correlated lattice problems as well as for 
non-equilibrium situations in mesoscopic systems. 

In section 2 we describe several exact properties of the auxiliary particle rep- 
resentation, while the conserving slave boson theory is developed and evaluated 
in section 3. 



2 Exact Auxiliary Particle Representation 

2.1 The SU(N)xSU(M) Anderson Impurity Model 

The auxiliary or slave boson method is a powerful tool to implement the 
effective restriction to the sector of Fock space with no double occupancy imposed 
by a large on-site repulsion U. The creation operator for an electron with spin a 
in the d-level is written in terms of fermionic operators f a and bosonic operators b 
as dl = fib. This representation is exact, if the constraint that the total number 
operator of auxiliary fermions f a and bosons b is equal to unity is obeyed. /J 
and b^ may be envisaged as creating the three allowed states of the impurity: 
singly occupied with spin a or empty. 

In view of the possibility of both FL and NFL behavior in quantum impurity 
systems mentioned in the introduction it is useful to introduce M degenerate 
channels for the conduction electron operators cj. , labeled fi = 1,2,..., M, in 
such a way that in the limit of impurity occupation number rid — * 1 (Kondo limit) 
the M-channel Kondo model is recovered, i.e. the model obeys an SU(M) channel 
symmetry. The slave bosons then form an SU(M) multiplet 6p which transforms 
according to the conjugate representation of SU(M), so that fx is a conserved 
quantum number. Generalizing, in addition, to arbitrary spin degeneracy N, 
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a = 1,2,..., TV, one obtains the SU(N)xSU(M) Anderson impurity model in 
slave boson representation 

where the local operator constraint Q = J2 a /I/u+E^ fysfyii = 1 must be fulfilled 
at all times. 



2.2 Gauge Symmetry and Exact Projection onto the Physical Fock 
Space 



The system described by the auxiliary particle Hamiltonian (2.1) is invariant 
under simultaneous, local U(l) gauge transformations, f a — > faC l ^^ T \ — > 
bp,e l ^ T \ with 0(t) an arbitrary, time dependent phase. While the gauge sym- 
metry guarantees the conservation of the local, integer charge Q, it does not 
single out any particular Q, like Q = 1. In order to effect the projection onto the 
Q = 1 sector of Fock space, one may use the following procedure j|, g]: Consider 
first the grand-canonical ensemble with respect to Q and the associated chemi- 
cal potential —A. The expectation value in the Q = 1 subspace of any physical 
operator A acting on the impurity states is then obtained as 

(A) = Urn 8C - L r = lim , (2.2) 

V 7 A^oo ^tr[e-/ 3 ( ff +^)] G A^oo(g) G ' 1 ' 

where the index G denotes the grand canonical ensemble and £ is the fugacity 
C = e _,3A . In the second equality of Eq. ( |2.2| ) we have used the fact that any 
physical operator A acting on the impurity is composed of the impurity electron 
operators d a , S a , and thus annihilates the states in the Q = sector, A\Q — 0) = 



0. It is obvious that the grand-canonical expectation value involved in Eq. (2.2) 
may be factorized into auxiliary particle propagators using Wick's theorem, thus 
allowing for the application of standard diagrammatic techniques. 

It is important to note that, in general, A plays the role of a time dependent 



gauge field. In Eq. (2.2) a time independent gauge for A has been chosen. In this 
way, the projection is only performed at one instant of time, explicitly exploit- 
ing the conservation of the local charge Q. Thus, choosing the time independent 
gauge means that in the subsequent development of the theory, the Q conser- 
vation must be implemented exactly. This is achieved in a systematic way by 
means of conserving approximations ||Io| , i.e. by deriving all self-energies and 
vertices by functional derivation from one common Luttinger-Ward functional 
$ of the fully renormalized Green's functions, 



E b ,f,c = 8*{G bt Gf,G B }/8G b ,f tC . 



(2.3) 
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This amounts to calculating all quantities of the theory in a self-consistent way, 
but has the great advantage that gauge field fluctuations need not be considered. 



2.3 Infrared Threshold Behavior of Auxilary Propagators 



The projection onto the physical subspace, Eq. ( |2.2| ), implies that the pseudo- 
fcrmion and slave boson Green's functions Gf,Gb are definied as the usual time- 
ordered, grand canonical expectation values of a pair of creation and annihilation 
operators, however evaluated in the limit A — > oo. It follows that the traces 
involved in Gf, Gb are taken purely over the the Q = sector of Fock space, 
and thus the backward-in-time contribution to the auxiliary particle propagators 
vanishes. Consequently, the auxiliary particle propagators are formally identical 
to the core hole propagators appearing in the well-known X-ray problem JTT[ , and 
the long-time behavior of Gf {Gb) is determined by the orthogonality catastrophe 
JL2| of the overlap of the Fermi sea without impurity (Q — 0) and the fully 
interacting conduction electron sea in the presence of a pseudofermion (slave 
boson) (Q = 1). It may be shown that the auxiliary particle spectral functions 
have threshold behavior with vanishing spectral weight at T — for energies 
u) below a threshold E Ql and power law behavior above E a , Af^ito) cx Q(lu — 
E Q )uj- a '- b . 

For the single-channel Anderson model, which is known to have a FL ground 
state, the threshold exponents may be deduced from an analysis in terms of 
scattering phase shifts, using the Friedel sum rule, since in the spin screened FL 
state the impurity acts as a pure potential scatterer [l^, |l6) , 

a/= 2n^n| ! « b = l-^ (N>l,M=l) (2.4) 

These results have been confirmed by numerical renormalization group (NRG) 
calculations |l7]] and by use of the Bethe ansatz solution in connection with 
boundary conformal field theory (CFT) [|l8). On the contrary, in the NFL case 
of the multi-channel Kondo model the threshold exponents have been deduced 
by a CFT solution jl|] as 

M N 

at = , a b = (N > 2, M > N) (2.5) 

' M + N ' 6 M + N V " ' " ' K ' 

Since the dependence of a/, ab on the impurity occupation number rid shown 
above originates from pure potential scattering, it is characteristic for the FL 
case. The auxiliary particle threshold exponents are, therefore, indicators for FL 
or NFL behavior in quantum impurity models of the Anderson type. 
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Figure 1 Diagrammatic representation of the Bethe-Salpeter equation defining 
the conduction electron-pseudofermion T-matrix t' c ^. 



3 Conserving Slave Particle T-Matrix Approximation 
3.1 Non-Crossing Approximation (NCA) 

The conserving formulation discussed in section 2.2 precludes mean field ap- 
proximations which break the U(l) gauge symmetry, like slave boson mean field 
theory. Although the latter can in some cases successfully describe the low T be- 
havior of models with a FL ground state, it leads to a spurious phase transition 
at finite T and, in particular, fails to describe NFL systems. 

Rather, the approximation should be generated from a Luttinger-Ward func- 
tional $. Using the hybridization V as a small parameter, one may generate 
successively more complex approximations. The lowest order conserving approx- 
imation generated in this way is the Non-crossing Approximation (NCA) |2(], [2l[] , 
defined by the first diagram in Fig. labeled "NCA" . The NCA is successful in 
describing Anderson type models at temperatures above and around the Kondo 



temperature Tk, and even reproduces the threshold exponents Eq. (2.5) for the 
NFL case of the Anderson impurity model. However, it fails to describe the FL 
regime at low temperatures. This may be traced back to the failure to capture 
the spin-screened Kondo singlet ground state of the model, since coherent spin 
flip scattering is not included in NCA, as seen below. 



3.2 Dominant Contributions at Low Energy 

In order to eliminate the shortcomings of the NCA mentioned above, we may 
use as a guiding principle to look for contributions to the vertex functions which 
renormalize the auxiliary particle threshold exponents to their correct values, 
since this is a necessary condition for the description of FL and NFL behavior, 
as discussed in section 2.3. As shown by power counting arguments p3], there 
are no corrections to the NCA exponents in any finite order of perturbation 
theory. Thus, any renormalization of the NCA exponents must be due to singu- 
larities arising from an infinite resummation of terms. In general, the existence 
of collective excitations leads to a singular behavior of the corresponding two- 
particle vertex function. In view of the tendency of Kondo systems to form a 
collective spin singlet state, we expect a singularity in the spin singlet channel 
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spin 
fluctuations 



charge 



Figure 2 Diagrammatic representation of the Luttinger-Ward functional gener- 
ating the CTMA. The terms with the conduction electron lines running clockwise 
(labelled "spin fluctuations" ) generate J 1 ^' , while the terms with the conduction 
electron lines running counter-clockwise (labelled "charge fluctuations" ) generate 
T^ cb \ The two- loop diagram is excluded, because it is not a skeleton. 



of the pseudofermion-conduction electron vertex function. It is then natural to 
perform a partial resummation of those contributions which, at each order in 
the hybridization V, contain the maximum number of spin flip processes. This 
amounts to calculating the conduction electron-pseudofermion vertex function 
in the "ladder" or T-matrix approximation, where the irreducible vertex 

is given by V 2 G b . The Bethe-Salpeter equation for T (c f> reads (Fig. Q), 

- V 2 T ^ Gbfi {iuj n + iw£ - x (3.6) 

where a, r, a' ', r' represent spin indices and fj, a channel index. A similar integral 
equation holds for the charge fluctuation T-matrix it is obtained from T^ c ^ 

by interchanging f a <-> and c CTM <-> cj.„. Inserting NCA Green's functions for 
the intermediate state propagators of Eq. ( |3.6[ ), we find at low temperatures 
and in the Kondo regime (n^~ 0.7) a pole of T {c ^ in the singlet channel as a 
function of the center-of-mass (COM) frequency fl, at a frequency which scales 
with the Kondo temperature, SI = f2 c / ~ — Tr-. Similarly, the corresponding T- 
matrix T( cf) ) in the conduction electron-slave boson channel, evaluated within the 
analogous approximation, develops a pole at negative values of £1 in the empty 
orbital regime (nrf~0.3). In the mixed valence regime (n^ ~ 0.5) the poles in 
both and coexist. The appearance of poles in the two-particle vertex 

functions T^ c ^ and T^ cb \ which signals the formation of collective states, may 
be expected to influence the behavior of the system in a major way. 
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(c) (d) 



Figure 3 Set of skeleton diagram parts not contained in CTMA which cancel 
in the infrared limit to leading and subleading order in the external frequency, 
w->0. 

3.3 Self-consistent Formulation 

The approximation considered so far is not yet consistent: The spectral weight 
produced by the poles of T^ c ^ and at negative frequencies Q is strictly pro- 
hibited in the limit T — > by the threshold property of auxiliary particle vertex 
functions (compare section 2.3). However, recall that a minimum requirement on 
the approximation used is the conservation of gauge symmetry. This requirement 
is not met when the integral kernel of the T-matrix equation is approximated 
by the NCA result. Rather, the approximation should be generated from the 
Luttinger-Ward functional shown in Fig. ||. It is defined as the infinite series of 
all vacuum skeleton diagrams which consist of a single ring of auxiliary particle 
propagators, where each conduction electron line spans at most two hybridiza- 
tion vertices. The first diagram of the infinite series of CTMA terms corresponds 
to NCA. By functional differentiation with respect to the conduction electron 
Green's function and the pseudofermion or the slave boson propagator, respec- 
tively, the shown i> functional generates the ladder approximations T^ c ^\ 
for the total conduction electron-pseudofermion vertex function and for the total 
conduction electron-slave boson vertex function (Fig. [l]). The auxiliary particle 
self-energies are obtained in the conserving scheme as the functional derivatives 
of <& with respect to Gf or G&, respectively (Eq. (|2.3|)), and are, in turn, non- 
linear functionals of the full, renormalized auxiliary particle propagators. This 
defines a set of self-consistency equations, which we term conserving T-matrix 
approximation (CTMA). 

The CTMA is justified on formal grounds by a cancellation theorem for all 
diagrams not included in the generating functional $ As seen in Fig. |2|, 
<f> includes all contributions where a conduction electron line spans up to two 
hybridization vertices; contributions not included contain at least one conduc- 
tion electron arch spanning four vertices. These contributions cancel pairwise at 
arbitrary loop order to leading and subleading order in the external frequency, 
as illustrated in Fig. ||. 
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Figure 4 The fermion and boson threshold exponents a/, at, are shown for 
N = 2, M = 1 in dependence of the average impurity occupation ra<j. Solid lines: 
exact values, Eq. (2.4); Symbols with error bars: CTMA; dashed lines: NCA. 



3.4 Results 

We have solved the CTMA equations numerically for a wide range of impurity 
occupation numbers rid from the Kondo to the empty impurity regime both for 
the single-channel and for the two-channel Anderson model down to tempera- 
tures of the order of at least 10 _2 Tr-. 

The solution of the CTMA equations forces the T-matrices to have vanishing 
spectral weight at negative COM frequencies O. Indeed, the numerical evaluation 
shows that the poles of T^ c ^ and T( cb ) are shifted to £1 = by self-consistency, 
where they merge with the continuous spectral weight present for fl > 0, and 
thus renormalize the threshold exponents of the auxiliary spectral functions. For 
AT = 2, M = 1 the threshold exponents otf, <Xb extracted from the numerical 
solutions are shown in Fig. ^|. In the Kondo limit of the multi-channel case 
(N >2,M = 2, 4) the CTMA solutions are found not to alter the NCA values 
and reproduce the the correct threshold exponents, a/ = M/(M + N), at, = 
N/(M + N). 

The good agreement of the CTMA exponents with their exact values over the 
complete range of for the single-channel model and in the Kondo regime of the 
multi-channel model may be taken as evidence that the T-matrix approximation 
correctly describes both the FL and the non-FL regimes of the SU(N)xSU(M) 
Anderson model (N=2, M=l,2,4). The static spin susceptibility \ of the single- 
and of the two-channel Anderson model in the Kondo regime calculated within 
CTMA as the derivative of the magnetization M = fiWBfa/T - n fl) with respect 
to a magnetic field H is shown in Fig. |5|. Good quantitative agreement with exact 
solutions is found for N = 2, M = 1 (FL). For N = 2, M = 2 (NFL) CTMA 
correctly reproduces the exact p5| logarithmic temperature dependence below 
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Figure 5 Static susceptibility of the single-channel (N = 2, M = 1) and the 
two-channel (N = 2, M = 2) Anderson impurity model in the Kondo regime 
(Ed = — 0.8D, T = O.ID, Lande factor g = 2). In the single-channel case, the 
CTMA andNCA results are compared to the NRG result ( T = 0, same pa- 
rameters) §3| and to the Bethe ansatz for the Kondo model gj]. The CTMA 
susceptibility obeys scaling behavior in accordance with the exact results (not 
shown) . 

the Kondo scale Tg. In contrast, the NCA solution recovers the logarithmic 
behavior only far below Tk- 



4 Conclusion 

We have reviewed a novel technique to describe correlated quantum impurity 
systems with strong onsite repulsion, which is based on a conserving formula- 
tion of the auxiliary boson method. The conserving scheme allows to implement 
the conservation of the local charge Q without taking into account time depen- 
dent fluctuations of the gauge held A. By including the leading infrared singular 
contributions (spin flip and charge fluctuation processes), physical quantities, 
like the magnetic susceptibility, are correctly described both in the Fermi and 
in the non-Fermi liquid regime, over the complete temperature range, including 
the crossover to the correlated many-body state at the lowest temperatures. As 
a standard diagram technique this method has the potential to be applicable to 
problems of correlated systems on a lattice as well as to mesoscopic systems out 
of equilibrium via the Keldysh technique. 

We wish to thank S. Bocker, T.A. Costi, S. Kirchner, A. Rosch, A. Ruckenstein 
and Th. Schauerte for stimulating discussions. This work is supported by DFG 
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